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Abstract 

A brief review of recent research on soliton and black hole solutions of Einstein's equations 
with nonlinear field sources is presented and some open questions are pointed out. 
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1 Introduction 

During the past few years there has been a lot of interest in sohton (by which I mean non- 
singular finite energy stationary solution) and black hole solutions of Einstein's equations 
with nonlinear field sources. In this brief survey I shall outline the most interesting, 
in my view, results of these studiesQ. Since my intention is to emphasize the model- 
independent aspects of the problem I shall restrict the details of specific models to a 
minimum (with one important exception of the Einstein- Yang-Mills theory which deserves 
separate analysis). 

It is rather surprising that although solitons in various nonlinear field theories in fiat 
space were intensively investigated in the past, the analogous problem in general relativ- 
ity received little attention until recently. Presumably this lack of interest in Einstein's 
equations with nonlinear field sources was due to two widely accepted beliefs (apart from 
the psychological factor due to the formidable structure of equations). First, it has been 
thought that the gravitational interaction, being very weak, cannot change qualitatively 
the spectrum of soliton solutions of a theory in which gravity was neglected. Second, 
according to the so called "no-hair" conjecture, a stationary black hole is uniquely de- 
termined by global charges given by surface integrals at spatial infinity, such as mass, 
angular momentum, and electric (magnetic) charge. It was believed that this property 
of black holes, following from the uniqueness [2] and the non-existence [3] results proven 
rigorously for various linear field sources, persists for general matter sources. 

The situation has changed radically when it was discovered that the Einstein- Yang- 
Mills (EYM) equations admit static non-abelian soliton [4] and black hole [5] solutions. 
This came as a surprise because i) neither the vacuum Einstein equations [6], nor pure 
Yang-Mills equations [7] have soliton solutions, and ii) "colored" black holes have non- 
abelian hair which is not associated with any conserved charge. These unexpected results 

have put into question the two beliefs mentioned above and have launched intensive 
tried to make this survey a complementary up-date of an excellent review on this subject written 
two years ago by Gibbons [1]. 
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investigations of Einstein's equations with nonlinear field sources. 

There are several, more or less independent, motivations to study this subject. Let 
us mention some of them. First and most important, this research may have serious 
conceptual implications. Our understanding of general relativity is to large extent based 
on the knowledge of few exact solutions {e.g. Kerr-Newman, Priedmann). However, as the 
results of [4] and [5] showed, the intuitions based on solutions with linear field sources 
fail in more general situations. In this sense the analysis of Einstein's equations with 
nonlinear field sources may shed new light on generic properties of solutions of Einstein's 
equations. 

Particularly interesting are non-perturbative effects of gravity, as, for instance, the 
existence of gravitational equilibria of non-abelian gauge fields [4] with finite energy. They 
result from the cancellation of gauge and gravitational singularities which implements the 
old idea due to Einstein that gravity may regularize ultraviolet divergencies in field theory. 
Such non-perturbative phenomena may have deep consequences at the quantum level. 

Second, although we do not have yet a unified quantum theory of all interactions, 
we obviously can couple gravity to the standard model at the classical level. Then, it is 
natural to expect that the effects of gravity become important when the relevant energy 
scale is close to the Planck scale, i.e. in the very early Universe. This fact is notoriously 
ignored in cosmological literature. Por example, in discussions how inflation solves the 
monopole problem the gravitational properties of monopoles are neglected, which might 
be a serious omission. 

Third, the long time behaviour of perturbed solitons is closely related to the main 
unsolved problem in general relativity: cosmic censorship hypothesis. 

Pinally, there is a purely mathematical motivation to study these problems. In the 
spherically symmetric case (which is physically most interesting) Einstein's equations 
with nonlinear field sources reduce to dynamical systems having very rich structure. It 
seems that modern methods of bifurcation and critical point theories may be succesfuUy 
applied to these systems. 

The rest of this brief review is organized into two sections dealing with solitons and 
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black holes, respectively. In Section 2 the general properties of gravitating solitons are 
described. For weak coupling the perturbative effects of gravity are examined, while for 
strong coupling the critical dependence on the coupling constant is analyzed. Next, the 
non-perturbative phenomena are discussed, in particular the Bartnik-Mckinnon solutions 
of the EYM equations. Section 3 is devoted to black holes with nonlinear hair and the 
current status of "no-hair" conjecture. This review is concerned with fundamental aspects 
of gravitating solitons and hairy black holes and I shall not discuss here an interesting 
problem of possible astrophysical and cosmological relevance of these solutions. 

I shall use units in which c = 1; all other dimensional parameters are explicitely 
written. 

2 Globally regular solutions 

Let us consider a nonlinear field theory in fiat space and suppose that the field equations 
have a soliton solution. For the purpose of this review I define the term soliton to mean a 
stationary solution which is everywhere non-singular and has finite energy^. The following 
question arises: What happens when the gravitational interaction is included into the 
model? In particular, how does gravity affect the spectrum of soliton solutions? To study 
this problem we have to consider Einstein's equations with a solitonic field as the source. 

Before doing this, it is worth recalling that a theory which satisfies the dominant 
energy condition and is scale invariant, cannot have soliton solutions. The reason is that 
under the scaling transformation x Xx, the energy scales a.s E ^ XE, so a stationary 
solution, being the extremum of energy, must have E = 0. Since the energy is positive 
this implies that the solution is trivial. Thus, our fiat space theory, having a soliton, 
must have a scale of length, call it Lq, which determines the characteristic size of the 
soliton (of course there is also a corresponding scale of energy). The gravitational cou- 
pling introduces into the model the additional dimensional parameter, Newton's constant 
G, which allows to define the second scale of length Lg (and energy). Thus typically the 
^ I do not require stability, so this meaning is broader than that usually used in literature. 
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Einstein-solitonic-matter system has two scales of length and the ratio a = Lg/Lg is 
a dimensionless parameter characterizing the model. Although my considerations are 
supposed to be model- independent, it is helpful to keep in mind a specific model for 
illustration. For this purpose I shall use the Einstein-Skyrme model [8,9]. In this model 
the two scales of length are given hy Lq — 1/ev and Lg — G^^'^/e, where v and e are 
dimensional parameters, and a — G^l'^v (cf.[8]). 



Weak coupling 



For small a one may apply the standard perturbative argument to prove the existence 
of gravitating solitons. The field equations may be expressed schematically in the form 

F(/,a) = 0, (1) 

where F is the differential operator and / denotes the collection of all unknown metric 
and matter field variables. Suppose that the parameters which determine Lq are fixed, 
hence the limit a — corresponds to switching-off gravity. By assumption, for a = 
(without gravity) there exists a soliton solution /q satisfying 

i^(/o,0) = 0. (2) 

A natural idea to show the existence of solutions of Eq.(l) is to use the implicit function 
theorem . In the present context the appropriate version of the implicit function theorem 
is the following. Let X, Y be Banach spaces, F : X x i?^ ^ y be a smooth mapping and 
let a point p = (/o, 0) e X x i?^ be a solution of Eq.(2). Then, if the derivative mapping 
DfF(p) is bijective, the theorem guarantees that in the neighbourhood of p there exists 
a solution f{a) of Eq.(l), such that /(O) = /q. The derivative DfF is just the linearized 
operator |j one obtains when linearizing F. In the most interesting physical situation 
the soliton /o is linearly stable which means that the eigenvalue problem 

SF 

jjip)C = ^'C (3) 
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has only positive eigenvalues, so the linear operator jj{p) is invertible, and, after choos- 
ing appropriate Banach spaces, bijective. If the spectrum of Eq.(3) is not positive (as 
it happens for unstable solitons), we have to make sure in addition that there are no 
zero modesQ. We conclude therefore that, in general, the soliton solution persists for 
sufficiently small a. Since the gravitational binding energy is negative, the energy of a 
gravitating soliton decreases with a. In "everyday" situations a is extremely small, so 
the gravitational effects are negligible. 

It might seem from the above argument that in the region of weak coupling nothing 
spectacular can happen. However, sometimes even for weak coupling interesting things 
occur. To see this, notice that I assumed above that a = corresponds to = (G = 0). 
However, since a = Lg/Lo, the limit a ^ may be obtained in another way, namely 
by keeping Lg fixed and taking Lq — > oo. This results in Einstein's equations with some 
truncated matter sources. For example in the Einstein-Skyrme theory in this limit {v — > 
with G fixed) the sigma-model term disappears from the lagrangian. An interesting 
situation arises when such a hmiting theory, albeit usually unphysical, has a soliton 
solution. Then, in analogy to Eq.(l) we have 

F(/,a) = 0, (la) 

and for a = there is a solution /o satisfying 

F(/o,0) = 0. (2a) 

The bar over F and / means that although Eqs.(l) and (la) are equivalent when a > 0, 
the operators F and F are different; Eq.(la) is obtained from Eq.(l) by a scaling trans- 
formation (depending on a). The Eqs.(2) and (2a) describe different limiting theories; the 
procedures of rescaling and taking the limit a — do not commute. If the operator fj(p) 
is bijective, the implicit function theorem guarantees the existence of a solution /(a) for 
■^In gauge invariant theories there are pure gauge zero modes which one has to remove before applying 
the imphcit function argument. This can be achieved either by fixing the gauge or by working with the 
space of gauge orbits. 
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sufficiently small a. Then for small a there are two distinct solutions, / and /, which 
are perturbations of /o and /o, respectively. It should be emphasized that although the 
solution / is obtained by the perturbative argument, its occurence, from the point of view 
of the original flat space model, is non-perturbative. This interesting phenomenon was 
observed and described in more detail in the Einstein-Skyrme model [8]. It seems that 
this effect is characteristic for a certain class of models possessing two scales of length. 

Strong coupling 

I have argued above that, as long as the dimensionless coupling constant a is suffi- 
ciently small, a flat-space soliton survives the gravitational coupling. Will such a gravitat- 
ing solution persist for large values of a? The gravitational interaction becomes important 
when a ~ 1 since then the size of the soliton is of the order of its Schwarzschild radius. 
For example the typical size of the skyrmion is R ^ 1/ev while its mass is M ~ f/e, so 
GM/R ~ a^. Hence, one might expect that for a ~ 1 the gravitating soliton becomes 
unstable with respect to the gravitational collapse and there is a critical value ~ 1 
beyond which no solitons exist. In fact, this expectation was con&med numerically in 
several models [8-12]. For example in the Einstein-Skyrme model ao ~ 0.2. It is plausible 
that this behaviour is generic for stable strongly gravitating soliton^. 

To understand better the existence of a critical value cto let us recall the standard tech- 
nique of proving existence of solitons in fiat space. Typically one starts from a topological 
argument showing that the configuration space has a nontrivial topological structure and 
splits into homotopy classes labeled by a winding number (topological charge). Then one 
shows that within a given homotopy class the energy functional is bounded from below 
by a positive constant proportional to the topological charge and therefore there exists a 
positive lower bound for energy 

Eo = inf E[f]. (4) 

''Unstable gravitating solitons may exist for all values of a [13]. 
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The idea is to prove that this bound is attained by a regular function /o which is thereby 
a stationary solution. Technically this is a three step procedure: 

1. Construct a minimizing sequence {/„}, where /„ belong to the space M of regular 
finite energy functions, such that limn^ao E[fn\ = Eq. The existence of such a 
sequence is guaranted by Eq.(4). 

2. Prove that the sequence {/„} has a limit /o G M. 

3. Finally, show that E[fo\ — hm„_»oo -Sf/n]. This step is nontrivial because usually 
E[f] is not continuous in M. 

The proofs of existence of topological sohtons along these hues were given for several 
models involving non-abelian gauge fields [14]. 

Now, let us return to gravitating solitons and follow the above procedure. First notice 
that as long as we deal with asymptotically flat solutions with topology ont — const 
slices, the topological arguments are still vahd for gravitating solitons (as we shall see 
below in the black hole context the situation may be different). Similarly, the mass is 
bounded from below by a positive constant. However, when one tries to repeat the steps 
1-3 in the case of nonzero a, it turns out that for large a a minimizing sequence {fn} 
has no limit in M. In other words, the minimizing solution /o is singular [11]. This shows 
that one should be careful with topological arguments. The topological argument is only 
a hint for existence of a solution and has to supplemented by the real proof, which, as 
the above example teaches us, might not be a "mere technicality". 

The behaviour of gravitating solitons in the strong coupling region is not yet well 
understood. In particular, a general argument for the existence of a critical value ccq 
is lacking. One might try to approach the problem along the lines of [15], where some 
universal upper bounds for the binding energy of static configurations were derived. In 
specific spherically symmetric models another approach seems to be easier. Impose regular 
initial data at r = and show that, when a is large, a local solutions cannot be extended 
to a global solution. Actually, this is exactly how such systems are studied numerically. 



8 



A closely related issue is the question of stability of gravitating solitons. Let us recall 
that the the static solution / is said to be linearly stable if all linear perturbations 
Sf{t) around it remain bounded (in a suitable norm) in time. The usual procedure of 
investigating linear stability is to assume that 6f = e*"^*^ , where ^ is time-independent 
and linearize time-dependent version of Eq.(l) to get 

()F 

77^ = oo'W^ , (5) 

where W is some positive weight function. The solution / is unstable if there exists at 
least one eigenmode ^ (satisfying appropriate boundary conditions) with negative eigen- 
value u"^, since then Sf will grow exponentially in time. Otherwise the solution is linearly 
stable. If the flat space soliton is stable then its gravitating counterpart will also be stable 
for small a, because the eigenvalues u;^ change continuously with a. The mixed numerical 
and analytical analysis shows that the eigenvalues a;^ decrease with a and for a —>■ ao the 
lowest eigenvalue u"^ tends to zero. This is not a coincidence, but follows from the general 
relationship between zero modes of the linearized operator and bifurcation points^. This 
connection is one of the basic results of the catastrophe theory which allows to study 
stability problems in the essentially non-dynamical way [16] . For gravitating solitons this 
fact was first observed in the Einstein-Skyrme model [8], where at ao the branch of fun- 
damental gravitating skyrmions /(a) merges with another branch of unstable solutions 
/(a) (discussed in the text below Eq.(lb)). 

Gravitational desingularization 

Above I have discussed Einstein's equations with solitonic matter sources. Now, con- 
sider a different situation and suppose that a field source has no soliton solution in fiat 
space. Can the gravitational coupling help in this respect, i.e. might there exist globally 
regular solutions in the coupled Einstein- non-solitonic- matter theory? First, notice that 
^ In the presence of symmetries there might be zero eigenvalues which have nothing to do with 
biforcations (Goldstone modes). 
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if the field theory is scale invariant then, according to the remarks above, the necessary 
condition for existence of a soliton in Einstein-matter system is that gravity breaks scale 
invariance. For this reason the Einstein-massless scalar or the Einstein-Maxwell equations 
cannot have a sohton, while the Einstein- Yang-Mills (EYM) equations might have, be- 
cause the EYM theory has a scale of length given by -s/G/e, where e is the gauge couphng 
constant. 

Even when a field theory has no globally regular solutions, it usually will have solutions 
which behave correctly at large distances (but are singular at the origin). A typical 
example is the monopole solution (electric or magnetic) in electrodynamics. It should be 
stressed that the leading asymptotic behaviour at spatial infinity of such solutions will not 
be altered by the gravitational coupling because gravity decouples from sources at infinity. 
However, gravity will modify the short distance behaviour of solutions, and in particular 
may regulate their singularities. This idea is traced back to Einstein who believed that 
ultraviolet divergencies in field theory will be eliminated in some unified nonlinear theory. 
Let us illustrate such a desingularization phenomenon with two examples. 

Consider a spherical shell of radius r with uniform charge and mass density. Let e be 
its total charge and mo its bare mass. Ignoring gravity, the energy is given by 

m(r) = mo H (6) 

r 

which, of course, diverges as r tends to zero. Taking into account the gravitational binding 
energy, this formula is replaced by 

m(r) = mo H . (7) 

r r 

Solving this equation for m(r) and taking the limit r ^ 0, yields 

mir^O) = ^. (8) 



This heuristic argument can be made rigorous by solving the constraint equations in 
the Einstein-Maxwell theory [17]. The result (8) may be viewed as non-perturbative 
cancellation of two infinite self-energies: positive electrostatic and negative gravitational. 
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To show another example of similar type let us consider a very simple theory in which 
gravity is modelled by a massless scalar field called dilaton (scalar gravity). Dilaton 
couples in a universal way to matter (with lagrangian L^) through the term e~^"''^L^, 
where a is the dilaton coupling constant. If we take the electromagnetic field as matter, 
the static Maxwell-dilaton equations 

d(e-'«'^ * F) = 0, V V + ^e-2«^F2 = (9) 

have a simple spherical solution with magnetic charge 1/e [18] 

eF = AsmM(p, a0 = ln(l + — ). (10) 

er 

Without a dilaton the energy density of the monopole diverges at r = as Tqq ~ 
whereas in the present case Too ~ e~'^°''^F'^ ~ 1/?^^, hence the total energy is finite (equal 
to 1/ea). Although the total energy of the dilatonic monopole is finite, the solution (10) is 
still singular at r = 0. Below I shall discuss the EYM theory where gravitational couphng 
leads to non-perturbative globally regular solutions. 



Bartnik-Mckinnon solutions 



Probably a single most important result in the study of gravitating solitons was the 
discovery of solitons in the EYM system by Bartnik and Mckinnon (BM) in 1988 [4]. 
Since then many properties of BM solutions were understood but still some questions 
remain to be answered. Below I discuss briefly the main features of BM solitons and 
point out open problems. 

The Einstein- YM coupled system is described by the action 

1 



G 



(11) 



where i? is a scalar curvature and T — dA + eA A A is the Yang- Mills curvature of 
a YM connection A which takes values in the Lie algebra of the gauge group Q. BM 
considered Q — SU{2). There are two dimensional parameters in the theory: Newton's 
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constant G (of dimension length/ mass) and the gauge coupling constant e (of dimension 

mass~^/'^length^^^'^). 

Let us consider static spherically symmetric configurations. It is convenient to 
parametrize the metric in the following way 



(is2 ^ -A^Ndt^ + N-^dr'^ + r'^{d^'^ + sin^ M^"^) , 



(12) 



where A and are functions of r. 

Assuming that the electric part of the YM field vanishes (actually this is not a re- 
striction because one can show, [19], that there are no globally regular static solutions 
with nonzero electric field) the purely magnetic static spherically symmetric SU (2)- YM 
connection can be written, in the Abelian gauge, as [20] 



eA = WTid'd + (cot 'Ot^ + WT2) sm'ddip 



(13) 



where {i = 1, 2, 3) are Pauli matrices and w is a function of r. The corresponding YM 
curvature is given by 



eJF = w'ridr A d-d + w'T2dr A siwddip — (1 — w'^)T2,d'd A sin-t^d^j , 



(14) 



where prime denotes derivative with respect to r. 

Inserting these ansatze into the action (11) yields the reduced lagrangian]^ (where 
S=IQtiS Ldt) 



_rA{l-N)+AU 



dr 



where 



U 



[l-w 
2^ 



2\2 



Variation of L with respect of if , A, and A^ yields the field equations!] 



{ANw')' + -^Aw{l - w^) = , 



(15) 



(16) 



(17) 



^ Strictly speaking, to obtain (15) one has to add a surface term to the action (11). 

^The principle of symmetric criticality [21] guarantees that the variation of S within the spherically 

symmetric ansatz gives the correct equations of motion. 
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[r{l - N)^ = 2GU , (18) 

A' = ^Aw'^ . (19) 

Note that, using Eq.(19), A may be eliminated from Eq.(17). 

These equations have three exphcit abehan solutions. The first two are the vacua 
w = ±1, A = N = 1 with zero mass (or Schwarzschild A = 1, N = 1 — 2GM/r). The 
third solution 

^ = 0, A^l, + (20) 

describes the Reissner-Nordstrom black hole with mass M and magnetic charge 1/e. 

In order to construct non-abelian solutions which are globally regular one has to 
impose the boundary conditions which ensure regularity at r = and asymptotic flatness. 
The asymptotic solutions of Eqs. (17,18) satisfying these requirements are 

±w=l-br^ + 0(/) , N =1- 4b'^r^ + 0{r^) (21) 

near r — 0, and 

±w^l-- + Oi\), 7V=1-- + 0(1) (22) 

near r = oo. Here 6, c, d are arbitrary constants. 

Notice that, for the asymptotic behaviour (22), the radial magnetic curvature, Br — 
T2,{1 — w'^) / , falls-off as 1/r^, and therefore all globally regular solutions have zero YM 
magnetic charge. 

Bartnik and Mckinnon gave strong numerical evidence that there exists a countable 
sequence of initial values hn in G N) determining globally regular solutions w„, yl„, A^„. 
The index n labels the number of nodes of the function w„. The sequence of masses M„ 
of these solutions increases with n and tends to = 1 for n — > oo (the unit of mass is 
given by 1/G^/^e). If one defines the "fine structure" constant 7^ — he^, then the mass of 
BM solitons is of the order mpi/j while their effective size is of the order lpi/'~f. Notice 
that the classical treatment of such objects makes sense physically only if 7 is small since 
then the Compton radius is much less than the effective radius of BM solitons. 

The BM discovery raised a number of questions. Let us enumerate the most important 
ones: 
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• Existence: The first rigorous proof of existence of BM solitons was given by SmoUer 
et al. [22]. Another proof was recently presented by Breitenlohner et al. [23]. Both 
proofs use the methods of dynamical system theory to analyze the a priori be- 
haviour of orbits starting with regular initial data at r = and find that there 
exists a countable family of connecting orbits (separatrices) which correspond to 
the BM solutions. Both proofs leave open the question of uniqueness of BM solu- 
tions (amongst globally regular solutions of Eqs.(17-19)). 

• Stability: The linear stabihty analysis of BM solutions was carried out by 

Straumann and Zhou [24]. They derived an eigenmode equation, such as Eq.(5), 
which governs the time evolution of small perturbations around the BM solitons and 
demonstrated that it has at least one negative eigenvalue (unstable mode). Later 
[25] they also analyzed the long time behaviour of perturbed BM sohtons and ar- 
gued that, depending on an initial perturbation, they either disperse to infinity or 
collapse to form the Schwarzchild black hole. 

• Raison d'etre: In the study of nonlinear equations perhaps even more important 

than proving existence of a solution is to understand the reason for its existence. 
The existence proofs mentioned above do not really explain what are the essen- 
tial features of the SU (2)-YM gauge group which are responsible for the existence 
of BM solutions and their basic properties, such as discreteness and instability. 
Soon after the discovery of BM solutions P. Mazur suggested (private commu- 
nication, 1990) that they should be regarded as gravitational sphalerons. Let us 
recall that sphaleron is a saddle point solution whose existence follows from the 
Lusternik-Snirelman mini-max (mountain pass) construction applied to the space 
of non-contractible loops (or paths joining topologically non-equivalent vacua) of 
field configurations. Sphalerons were first discovered by Taubes and Manton in 
spontanously broken gauge theories (YM-Higgs) [26], but they exist in many other 
models which possess non-contractible loops in configuration space [27]. The first 
published attempt [28] to interpret the BM solutions as sphalerons was unsuccesful. 



because the authors considered loops with too weak decay condition for the YM 
connection which led to a disaster of having the infimum of energy over all loops 
equal to zero (c/. footnote 9). 

Recently Sudarsky and Wald (SW) [29] proposed, in the language of Hamilto- 
nian formulation of general relativity, a very interesting modification of the original 
Taubes-Manton argument. In the case of n = 1 BM solution, the SW argument 
is, in essence, equivalent to the mini-max procedure for paths joining two topolog- 
ically inequivalent vacua. However, in contrast to the mini-max construction, the 
SW mechanism can be naturally extended to account also for the existence of n > 1 
BM solutions. The SW argument got support by the discovery of solutions similar 
to BM solutions in another theory involving SU{2)-YM field [18]. It is very unlikely 
that the SW argument it its present form can be converted into a genuine proof. 
However, in the spherically symmetric case the situation simplifies considerably and 
seems to be tractable rigorously. I shall come back to this point when discussing 
the YM field on Schwarzschild background. 

Spectrum of masses: The discrete spectrum of masses of BM solutions has an 
accumulation point (upper bound) at M = l/(eG^/^). This follows from the fact 
that for n — > oo the BM solutions tend (non-uniformly) to the extremal Rcissncr- 
Nordstrom solution with magnetic charge equal to 1/e. Can one construct an in- 
equality M < l/{eG^^'^) which is saturated by the hmiting solution? 

Spherical symmetry: In the EYM system the Birkhoff theorem is not valid since 
a spherically symmetric configuration need not be static. Is the converse result true, 
i.e. does staticity imply spherical symmetry? Showing this would be the first step 
towards the proof of uniqueness of BM solutions (of static SU {2)-YM equations). 
One could try to attack this problem using the technique of Bunting and Masood-ul- 
Alam [30] of finding a conformal transformation of spatial 3-metric into a zero-mass 
metric with non-negative scalar curvature. As far as I am aware (W. Simon, private 
communication) this approach encounters serious technical difficulties. A closely 
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related question: are there stationary axisymmetric analogues of BM solutions? 

• Desinguleirization: The BM solutions are non-perturbative in the gravitational 
constant G - nota bene in the EYM theory there is no dimensionless parameter 
which might serve as a perturbative parameter (changing the coupling constants e 
and G only changes the scale). Although the flat-space static YM equations have 
no globally regular solutions [7] , they possess spherically symmetric solutions which 
at large distances behave similarly to BM solutions. These are solutions of Eq.(17) 
with A = N = 1 [31]. Their asymptotic behaviour at infinity is given by Eq.(22), 
whereas near r = they behave as 

w ~ -\A'sin(-— Inr + const) . (23) 

so the YM curvature T is singular at r = and the total energy is infinite. I beheve 
that the BM solutions may be viewed as the regularized (by gravity) version of these 
singular flat-space solutions. In a sense the proof of Breitenlohner et al. implements 
this idea, because they start from the phase portrait of Eq.(17) in flat space and 
analyze how it is modified by gravity. It would be of interest to pursue this idea 
further. One possible direction is to search for exact form of BM solutions, since 
having it wc probably could sec how the scale of length G^^'^/e provides a cut-off 
for singular behaviour (23). Recent results in the EYM-dilaton theory (Bizoh and 
Chmaj, work in progress) suggest that the hope of finding an analytical expression 
for BM solutions is not unrealistic. 

3 Black holes 

Let us now turn to another interesting class of asymptotically fiat solutions, namely black 
holes. We restrict attention to the region outside the horizon. Then, in the stationary case, 
Einstein's equations reduce to the elliptic boundary value problem between horizon and 
infinity. 
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There are more black hole solutions than globally regular ones. This is due to two 
important differences between these two cases. First, stationary globally regular solutions 
are extrema of mass, while stationary black holes extremize the mass provided that the 
area of the horizon is fixed (and the angular momentum is fixed, if there is nonzero 
rotation)! [29]. This follows from the first law of black hole dynamics 



where M, k, A, J, and f2 are the black hole mass, surface gravity, area, angular momentum, 
and angular velocity of the horizon. In other words the presence of horizon introduces 
the boundary condition which breaks scale invariance. This allows for nontrivial black 
hole solutions in scale invariant theories, as for instance the Schwarzschild solution. 

Second, solutions which are singular in fiat space may give rise to perfectly regular hair 
on black holes if the singularities are hidden behind the horizon. The Reissner-Nordstrom 
black hole is a typical example. 

In analogy to the discussion of globally regular solutions we can distinguish several 
situations depending on scaling properties. The first possibility is that the model is scale 
invariant, and the scale is given solely by the area of the horizon. In this case we have a 
continuous family of black holes, such as Schwarzschild or Reissner-Nordstrom solutions. 

The second possibility is that the corresponding fiat space model has a soliton and 
a scale Lq. The coupling of gravity gives a second scale Lg and the presence of horizon 
provides a third scale th, where vh is the radius of the horizon (for simplicity I consider 
spherically symmetric case). Thus, in analogy to Eq.(l), we can write the field equations 
in the form 



where a = Lg/ Lq and (3 = th/Lq. Note that the status of parameters a and [3 is different - 

a is a combination of coupling constants while (3 is determined by the boundary condition 
^In both cases there might be additional contributions to the variation of mass coming from other 
extensive parameters (such as e.g. charge). Mass is extremized at stationary solutions if these parameters 
are fixed [29]. 



G5M = —k5A + QSJ 
Stt 



(24) 



F{f,a,P) = 



(25) 
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(in the nonspherical case it is even not clear how to define (3). Nevertheless, as long as 
Th > 0, by rewriting the equations in the variable x = t/th, (3 is eliminated from the 
boundary condition and appears as a second (apart from a) dimensionless parameter in 
the equations . It is very convenient to think about Eq.(25) as describing black holes 
when /3 > 0, and globally regular solutions when /3 = 0. That is, we require / to belong 
to a space X of asymptotically fiat configurations which for /3 > satisfy the black-hole 
boundary conditions at r/j, while for = (when black hole becomes singular since the 
horizon shrinks to zero) are regular at the origin. The advantage of working with X is 
that the fiat-space solution /o, satisfying F(/o, 0, 0) = 0, belongs to X and may be used 
as the basis of the perturbative argument. Such an argument was proposed by Kastor 
and Traschen [32] who argued, using the implicit function theorem, that for sufficiently 
small a and j3 in the neighbourhood of /o there exist black hole solutions of Eq.(25). 
Although their heuristic argument is probably basically correct, it certainly should be 
pursued by more precise analysis which would take care of some technical subtleties, 
such as the effects of degenerate ellipticity at th- Anyhow, this reasoning is strongly 
supported by numerical results in several models having soliton solutions, where black 
hole counterparts with small horizon (sometimes referred to as black holes inside solitons) 
were found [8-12,33]. When /? ^ these black holes tend to the corresponding globally 
regular solutions. 

It turns out that such black holes exist only in a bounded region of the (a, /3)-plane. 
As was shown before, for /3 = there is a maximal value beyond which no globally 
regular solutions exist. Numerical and analytical analyses show that when jS grows, the 
critical value q;o(/9) monotonically decreases and goes to zero for some /?o [8-12,33]. Thus 
there is an upper bound for the radius of horizon. The heuristic explanation is that when 
the horizon gets larger, more and more of solitonic hair is swallowed and finally, when the 
radius of the horizon is comparable to the size of the sohton (/3 ~ 1), the whole soliton 
disappears inside the black hole (no "cloud" outside the horizon is left). The technical 
reason of non-existence of black holes with solitonic hair for large /? is the following. 
The standard no-hair proof (see papers by Bekenstein [3]) consists in constructing an 
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identity which has the form of a divergence equahng a non-negative definite quantity. 
After integrating this identity between horizon and infinity the divergence term vanishes, 
hence the integral of non-negative quantity is zero which imphes that the field (hair) is 
trivial. For a nonlinear source one can also construct a similar divergence identity but 
usually its right side has both positive and negative contributions, hence the proof fails. 
However in some cases one can show that when /3 is large enough the positive terms take 
over so the right side of the divergence identity is non-negative and the no-hair proof 
goes through. 

Note that the hmiting case a — corresponds to the decouphng of gravity. Thus to 
find a maximal value /3o it is sufficient to study solitonic fields on fixed Schwarzschild 
background. For skyrmionic hair this was done in [34]. 

The third possibility is that the fiat space model has no solitons, and gravity brings 
a scale Lg. Then the field equations have the form F{f,j) = 0, where 7 = rn/Lg. This 
situation arises in the EYM theory so let us now consider black holes in this model. As 
above, it is convenient to treat the limit 7 = as corresponding to a regular solution. As 
I have already described, in this case there is a countable family of BM solutions; let us 
denote them by /° (superscript refers to 7 = 0). For sufficiently small 7 we can repeat 
Kastor and Traschen's heuristic argument to predict the existence of a countable family 
of black holes fn{l), which are perturbations of the corresponding BM solutions. These 
solutions, so called colored black holes (CBHs), were found numerically independently 
by several authors [5] soon after the BM discovery. Recently, the existence of CBHs was 
established rigourously [22,23]. Both numerical and analytical results show that there is 
no upper bound for 7, i.e. for every 7 there is a countable family of CBHs, fn{j)- It should 
be stressed that CBHs with different 7 are essentially different; in particular they are not 
related by any simple scaling (contrary to some statements scattered in literature) . When 
7 ^ 00 the YM and gravitational fields decouple and we end up with the YM field on 
fixed Schwarzschild background. This model is extremely simple, nevertheless it seems 
to contain all the essential features which are responsible for the existence of BM-type 
solutions. For this reason it is an ideal ground for converting heuristic arguments (such as 
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the SW argument) into rigorous results. To my knowledge this model was not described 
in literature so I take this opportunity to give some details. 

SU(2)-YM on Schwarzschild background 

Let us rewrite Eqs. (17-19) using a dimensionless coordinate x = r/rH and take the 
limit 7 = erujG^^'^ — ^ cxd. In this limit the right side of Eqs. (18,19) vanishes, hence they 
are solved by the Schwarzschild metric 

= 1 - i, A = \, (26) 

X 

and therefore Eq.(17) reduces to 

((1-^K)' + ^^1-^')-0. (27) 



This equation follows from the variation of the energy functional 



dx . (28) 



It was shown first numerically and recently analytically (Bizoh, unpublished) that Eq. (27) 
has a countable family of solutions {tWn} (^ £ -^) which are regular for x > 1. The first 
member of this family is the vacuum = 1, for which the energy functional attains the 
global minimum E\wq\ = 0. For a discrete set of initial values w„(l) G (0, 1) there exist 
essentially non-abelian solutions which oscillate (n— l)-times within a strip (—1, 1) and 
tend to (—1)" for x — > oo. The discrete spectrum of energies Eq — 0,Ei 0.4795, E2 ~ 
0.4994, £'3 0.4999,... has an accumulation point at E^o — 1/2. This follows from the 
fact that for n — > 00 the solutions Wn tend pointwise (on any finite interval [1, Xq)) to the 
abelian solution Wqo = with energy E[woo] = 1/2. Using Eq.(27) one can easily show 
that on-shell the energy is given by 

— {l-w^)dx, (29) 
which proves that the solution Woo — has maximal energy. 
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For completeness let us note that the solution wi is known analytically [35] 



c — X 



(30) 



Wi = 



x + 3{c-l) 



where c= (3 + V^)/2. 

Notice that our system has Z2 symmetry, hence along with {wn} there is a corre- 
sponding family of reflected solutions {— The limiting solution Woo is a fixed point 
of reflection. 

The Hessian of the energy functional is given by 



A very interesting feature (found numerically) is that the Hessian evaluated at a solution 
Wn has exactly n negative eigenvalues. This property is a strong hint for a mechanism 
which gives rise to the family {w„}. Before discussing it let me make the following clari- 
fying remark. 

Remark. The key feature of the SU{2)-YM gauge group is the existence of large 
gauge transformations, i.e. topologically inequivalent cross sections of the YM bundle. 
They arise as follows. Choose the gauge Aq = and consider the spatial YM potential 
Ai on the t = const slice S of topology. The pure gauge configurations have the form 
Ai = diUU~^, where f/ is a map from S to the group manifold SU{2). Demanding that 
U 1 (sufficiently fast0), this becomes a homotopically nontrivial map U : ^ S^, 
classified by integer winding numbers. Thus, there is a countable set of topologically 
inequivalent vacua which cannot be continuously deformed into one another. This fact 
is the basis for the Taubes-Manton mini-max construction of sphalerons and for the 
Sudarsky-Wald construction of BM solutions. However, in the black hole setting the 
situation is different. The point is that now the space S has topology S*^ x R^ and the 
mapping U : S ^ S*^ is topologically trivial (it is not possible to compactify S by 

demanding U to be identity on the horizon). 

^It is essential that rdiU — > (and therefore rAi — > 0) as r 00. Otherwise the YM bundle is 
topologically trivial, as one can easily demonstrate by constructing curves of connections (with Ai — 
0{l/r) decay) joining large gauge transformated vacua. I thank R. Bartnik for pointing this out to me. 



rco r 1 1 

6'F{wm,0 = (1 - + -{5w' - l)e dx 

Jl I X x^ 



(31) 
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Although the topology of the SU{2)-YM field on the Schwarzschild background is 
trivial, the presence of two distinct vacua and Z2 symmetry seems to be sufficient for 
making the mini-max construction. Let F be a space of functions w for which the energy 
functional (28) is finite. Consider all paths in F connecting wq and —Wq. Since Wq and 
—wq are global minima of energy, on each path there is a point of maximal energy. The 
infimum over these maximal energies gives a saddle point (with exactly one unstable 
mode), which we identify as wi. It is tempting to repeat this procedure using wi and its 
refiection —Wi. However, this fails because there are paths joining wi and —Wi which go 
below the energy level One could try to remedy this obstacle by defining a space 

Fi consisting of all maxima on paths joining wq and —Wq, and making mini-max on Fi. 
Unfortunately, it is very doubtful that Fi will be sufficiently well-behaved to allow such 
construction]^, in particular almost surely Fi won't be a surface in F. Recently R. Wald 
suggested (private communication) that to obtain higher n solutions one should apply 
mini-max method not for paths but for higher dimensional objects, such as n-spheres 
which are invariant under refiections. 

Another very interesting approach was pursued by O. Popp (unpublished). The idea 
is to show that the energy functional is a Morse function on F \ Woo, i.e. that all critical 
points of E[w] (except Woo) have finite index (in the language of Morse theory, index = 
number of negative eigenvalues of the Hessian). To prove this one approximates E[w] by 
some truncated functional -FAf""^], defined by certain integral over the interval [1, A]. For 
A close to one, F\[w] has only three critical points: ±wo = 1 and Woo = and all these 
points have zero index. When A increases new critical points appear. At each bifurcation 
point A„ the solution Woo "sheds" a pair of new solutions (w„,— w„), and at the same 
time it acquires one additional negative eigenvalue {i.e. the index of Woo increases by 
one). Near the bifurcation point A„ the solution Wn is close to Woo, and therefore it has 
the same index (equal to n). As X increases the index of Wn remains constant (this fact 
is not yet proven rigorously but there is strong evidence that it is true). Showing that 
thank R. Wald for pointing this out to me. 
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Fx[w] E[w] for A — ^ oo concludes Popp's argument. 
No-hair conjecture 

The "no-hair" conjecture belongs to the folklore of general relativity since the late 
sixties (for a list of references, see [36]). The conjecture concerns the possible forms of 
stationary black holes. The idea is that the only classical degrees of freedom of black 
holes are those corresponding to non-radiative multipole moments, because "everything 
that can be radiated away will be radiated away" (cf. [37]). For a massless bosonic field 
with spin s, radiative multipoles have moments / > s, hence for pure gravity {s = 2) 
the conjecture allows for a monopole (mass) and a dipole (angular momentum) degrees 
of freedom. If there is an electromagnetic field (s = 1), the electromagnetic monopoles 
(electric and magnetic charge) are also allowed. The conjecture exludes the massless scalar 
field (s = 0) and all massive fields, because for these fields all multipoles are radiative. 

At the time when the "no-hair" conjecture was formulated this picture was perfectly 
consistent with the fact that the only known stationary black hole was the Kerr-Newman 
solution and nobody doubted that this was a unique stationary black hole solution of the 
Einstein-Maxwell equations (although this fact was established rigorously later on [2]). 
Morever, the conjecture was supported by several no-go results which showed that such 
fields as massless and massive scalar, massive vector or spinor cannot reside on stationary 
black holes [3]. 

The "no-hair" conjecture was formulated rather vaguely and therefore admits many 

interpretations. Nowadays most people agree that it should be meant as a statement 

concerning the uniqueness of stationary black holes. Let us adopt this viewpoint and try 

to formulate the "no-hair" conjecture more precisely^ For this purpose let me define a 

term global charge to mean a conserved charge associated with a massless gauge field, 

which is given by surface integral at spatial infinity. Mass, angular momentum, and 
^^I am concerned here with what the "no-hah" conjecture reaUy means, and not with mathematical 
assumptions underlying the uniqueness theorems (cf. [38]). 
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electric and magnetic charges are examples of global charges. Then the strongest variant 
of "no-hair" conjecture is 

Version 1. A stationary black hole is uniquely described by global charges. 

For a long time no counterexample to this version was known. However, recent dis- 
coveries of various black holes with nonlinear matter fields made clear that one cannot 
speak about uniqueness without specifying the sources. To see this consider an example of 
Einstein-Maxwell-dilaton theory. In this model Gibbons found very interesting black hole 
solutions [39]. In the so called string- inspired case (where the dilaton coupling constant 
a is equal to VG), the electrically charged solution is given by 

F^e'^^%dtAdr (32) 

where is a dilaton. It was shown recently that this (and the analogous magnetically 
charged solution) is a unique static black hole solution in this model [40]. However, the 
dilaton field is not associated with any conserved charge, and therefore at infinity this 
solution is indistinguishable from the Reissner-Nordstrom solution with the same mass 
and charge. Thus an observer at infinity would not be able to determine a solution by 
measuring global charges, unless she knows the matter content of the world. This leads 
us to the following modification: 

Version 2. Within a given model a stationary black hole is uniquely determined by 
global charges. 

Although in many cases, such as the electrovacuum, this version is actually a rigorous 
theorem, there are theories in which it is false. For example, as I discussed above, in the 
EYM theory there is a countable family of colored black holes. They carry no charge 
so their non-abelian field leaves no imprint at infinity. For a given mass (which is the 
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only global parameter) there are infinitely many CBHs with different areas of the horizon 
which clearly violates Version 2. However, this is not a physically serious counterexample, 
because CBHs are unstable [41] - under small perturbation they loose non-abelian hair 
and decay into Schwarzschild (for such unstable hair the word wig seems to be more 
appropriate). We should therefore add to Version 2 the following qualification: 

Version 3. Within a given model a stable stationary black hole is uniquely determined 
by global charges. 

Of course in the EYM theory this unique stable solution is Schwarzschild. Unfor- 
tunately, there are models which violate even this version [8,33]. For example, in the 
Einstein-Skyrme model, for a given mass there are three black holes, two stable and one 
unstable. One of the stable black holes is Schwarzschild while another one has skyrmionic 
hair. Since the skyrmionic hair is not connected to any conserved charge (and morever is 
topologically trivial), the existence of these two stable black holes provides a counterex- 
ample to Version 3. 

There seems to be no further possibility of relaxing Version 3 to accomodate such 
cases as the Einstein-Skyrme model. We conclude therefore that there is no universally 
valid formulation of the "no-hair" conjecture. 
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